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Executive  Summary 


Objective 

The  capability  to  numerically  model  the  scattering  of  sound  from  a  body  in  a  dispersive 
environment  has  been  developed  and  an  evaluation  has  been  made  of  the  applicability  of 
some  simple  physical  approximations.  The  model  is  based  on  a  rigorous  coupling  of  a  finite 
element  formulation  of  target  dynamics  with  a  Helmholtz  integral  equation  formulation  of 
fluid  loading.  The  effect  of  neglecting  multiple  scattering  and  Fresnel  effects  is  determined. 

Results 

It  has  been  determined  that  in  cases  of  practical  interest  where  a  acoustic  scatterer 
is  small  relative  to  the  dimensions  of  a  duct  in  which  it  lies,  multiple  scattering  effects  are 
small  except  when  the  target  is  within  approximately  its  own  diameter  of  the  boundary. 
Furthermore,  the  neglect  of  Fresnel  terms  has  small  effects.  Neglecting  these  two  physical 
effects  has  results  in  an  approximate  formulation  of  the  problem  which  is  more  efficiently 
solved  by  computer. 

Recommendations 

None. 
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INTRODUCTION 


An  approach  to  the  problem  of  calculating  the  pressure  field  scattered  by  a  structure 
immersed  in  a  nonuniform  fluid  environment  has  been  reported  in  McDaid  et  al.  (1992). 
The  method  is  based  on  the  Helmholtz  integral  equation  formulation  of  fluid  loading  in  a 
stratified  medium  and  on  a  finite  element  formulation  of  the  dynamics  of  the  structure. 

A  set  of  Helmholtz  Integral  Equations  has  been  derived  which  have  as  kernels  the 
Green’s  function  (and  its  gradient)  for  a  refractive  environment  with  boundaries.  The  rel¬ 
evant  equations  are  reproduced  below.  A  detailed  derivation  of  the  equations  has  been 
relegated  to  the  appendix  to  preserve  the  continuity  of  the  main  body  of  the  report. 

First,  there  is  the  Surface  Integral  Equation  relating  the  distribution  of  pressures  and 
velocities  on  the  wetted  surface  of  the  target,  which  have  been  induced  by  the  incident 
(ensonifying)  pressure  generated  by  a  distant  point  source  of  sound: 

surface  incident 

pressure  pressure 

1  — '  ' - ' 

-  P{Rsg\Rsp)  =  G{Rsp\Rsg)  + 

f  [P{Rsh\Rsp)'^R,h  ^s{Rsh\Rsg) 

JSsh.  ' - V - '  “  ' - V - ' 

surface  surface 

pressure  Green’s  fen 

where  Rsp^  Rsg^  and  Rsh  are  the  positions  of  the  locations  of  the  ensonifying  source  point, 
the  field  point  on  the  target  at  which  pressures  are  to  be  calculated  and  an  integration  point 
(dummy  variable)  on  the  target  surface.  The  quantity  fish  is  the  unit  outward  normal  vector 
to  the  wetted  surface  of  the  target. 

The  foregoing  equation  (1)  employs  an  approximation  that  reduces  the  computational 
burden,  viz.,  that  the  kernel  of  the  integral  equation  is  an  approximate  propagation  Green’s 
function,  G5,  rather  than  the  required  refractive  Green’s  function  G,  The  reason  for  using 
this  approximation  is  that  a  tractable  implementation  of  an  algorithm  for  a  normal  mode 
representation  of  the  Green’s  function  in  the  near  field  is  currently  not  available.  McDaid 
et  al.  (1992)  detailed  the  results  of  a  study  in  which  the  free-space  Green’s  function  Gq 
was  used  for  G5.  The  justification  of  the  use  of  the  free-field  approximation  is  that  the 
distances  separating  pairs  of  points  on  the  surface  are  small  compared  to  the  distances  over 
which  significant  refractive  effects  due  to  a  variable  sound  speed  will  occur;  and  the  interpoint 
distances  are  small  in  comparison  to  the  distance  to  the  surface  and  hence  to  target  “image” 
points.  This  latter  assumption  requires  that  only  those  configurations  be  investigated  for 
which  the  target  is  not  “near”  the  boundaries  of  the  medium. 


—  Gs{Rsh\Rsg)  Rsh  P{Rsh\Rsp)]  *  flshdS  ,  (1) 

surface  surface 

Green’s  fen  pressure 
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MULTIPLE  SCATTERING 


The  most  important  of  the  effects  that  were  ignored  by  McDaid  et  al.  (1992)  was 
multiple  scattering.  This  phenomenon  will  be  strongest  when  the  target  is  very  near  one  of 
the  boundaries  of  the  fluid.  In  the  present  study,  a  Neumann  series  representation  of  the 
Green’s  functions  appearing  in  the  Surface  Integral  Equation  (1)  is  used  to  approximately 
account  for  multiple  scattering  from  the  pressure  release  upper  surface  and  from  the  lower 
boundary.  The  series  can  be  used  to  generate  a  sequence  of  successively  more  accurate 
approximations  to  the  true  solution.  The  sequence  of  approximate  solutions  is  developed  in 
the  following  analysis. 

The  equation  for  the  Green’s  function  can  be  rewritten  in  the  form 

G{Robs\Rsp)  —  Go{Robs\Rsp)^ 

/  [Go{RuPPER\Rsp)^RuppER^i^obs\^UPPER)]  '  RfldS 

•^Supper 

+  /  [Go{RLOWER\Rsp)'^RLowERG{Robs\RLOWER)  + 

^Slower 

G{Robs\RLOWERy^  RlowerGo{RlOWEr\R  sp  )]  •  fijids 

-  I  [(eM-kl)Go(R,\R.,)G{R,\K)\<‘V  ^ 

J  u 

where  is  substituted  for  the  term  R^g. 

The  foregoing  equation  for  the  Green’s  function  can  be  written  in  the  form 


G  =  Go  +  £  {Go;  G}  , 


where  the  operator  C  (bilinear  in  Go  and  G)  represents  the  residual  effect  of  boundary 
scattering  and  refraction. 

The  problem  of  solving  for  the  Green’s  function  is  hence  cast  into  the  form  of  solving 
the  foregoing  fixed-point  problem.  Under  certain  circumstances  (i.e.,  when  the  right-hand 
side  is  a  contraction  mapping  and  a  starting  point  can  be  found  in  the  domain  of  attraction), 
the  problem  can  be  solved  by  the  method  of  Picard  iteration;  hence,  one  generates  a  sequence 
of  approximations,  G[n],  to  the  function  G  with  the  hope  that  the  sequence  converges  to  the 
desired  function.  Specifically,  one  starts  with  Go  and  one  constructs  the  sequence 
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G[o]  =  Go 

G[i]  =  Go->c{Go;G[o]} 
G[2]  =  Go  —  G  |Go;  G[i]| 


[n+l] 


Go-G{Go;G[„]} 


This  sequence  of  approximations  to  the  Green’s  function  will,  under  certain  circum¬ 
stances  (unknown  at  this  point),  converge.  In  this  context,  the  baseline  approach  to  be  used 
in  the  current  study  can  be  thought  of  as  the  0-th  order  approximation,  i.e.,  the  approx¬ 
imation  G[o]  is  used  in  the  kernel  of  equation  (1).  While  the  higher  order  approximations 
have  not  been  used  directly  in  this  study,  the  value  of  this  formulation  is  that  it  provides  a 
constructive  approach  for  calculating  each  term  of  the  sequence.  In  addition,  this  formula¬ 
tion  of  the  problem  provides  an  analytical  framework  for  the  expressions  used  to  generate 
estimates  of  multiple  scattering  effects. 

A  variant  of  the  lowest  order  approximation,  which  is  of  particular  interest  in  the  case 
of  a  duct  with  a  hard  bottom,  is 


Gs  =  Gio^]{Robs\Rsp)  =  Go  +  Go+  -  Go-  (2) 

This  approximation  will  capture  the  dominant  multiple  scattering  effects  for  a  target  near  the 
pressure  release  upper  boundary.  The  term  Go+  is  the  first  image  with  respect  to  the  bottom, 
and  Go-  is  the  first  (out  of  phase)  image  with  respect  to  the  top.  This  expression  represents 
a  special  subsequence,  which  has  been  summed  in  the  Neumann  series.  The  refractive  terms 
and  the  more  distant  images  (an  infinite  number  exist)  have  been  ignored. 

In  addition  to  the  equation  for  fluid  loading,  an  equation  relating  the  surface  pressures 
and  velocities  of  the  target  (possibly  elastic)  is  needed.  Although  a  shell  approximation  for 
the  linear  elastic  equations  describing  the  target  has  been  implemented,  only  rigid  targets 
are  used  in  the  present  study.  The  relevant  equations  have  been  described  by  Schenck  and 
Benthien  (1989)  and  are  not  reproduced  here.  Instead,  this  result  is  simply  expressed  as  the 
following  operator  relationship  between  the  functions  describing  the  pressure  and  velocity: 

Vp  •  Hsh  =  ^(p)  on  S  , 


where  S  is  the  wetted  surface  of  the  target.  The  foregoing  equations  are  solved  simultaneously 
to  determine  the  distribution  of  velocities  and  pressures  on  the  surface.  In  the  present 
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instance,  the  equations  reduce  to  the  condition  that  the  velocities  normal  to  the  target  are 
zero. 


Once  the  distributions  of  surface  pressures  and  velocities  have  been  calculated,  the 
scattered  field  is  computed  using  the  following  equation 


complete 

pressure 

field 

P{Robs\Rsp)  - 


source 

field 

contribution 

^  s 

G{Robs\Rsp)  + 


refractive  refractive 
surface  Green’s  Green’s 

pressure  function  function 


surface 

pressure 


target  field  contribution 


where  Robs  is  the  location  in  the  field  at  which  the  pressure  is  to  be  calculated. 


The  propagation  Green’s  function  is  given  a  “normal  mode”  representation  as 

*  OO  /  _ 

C(r/,  Zf\r,,  2s)  =  -  [r\lkl  - \A  (j)j{zf)(f)j{zs)  ,  (3) 

^  i=i  ^ 

where  (vs^Zs)  and  are  the  cylindrical  coordinates  (range  and  depth)  of  the  source 

point  and  field  point,  respectively,  of  the  Green’s  function. 

This  equation  (3)  is  substituted  into  the  equations  for  the  scattered  field  and  modal 
scattering  amplitudes  are  computed.  The  complete  scattered  field  for  a  target  can  be  com¬ 
puted  as  a  weighted  sum  of  modes.  In  those  cases  for  which  the  Green’s  function  cannot  be 
represented  as  the  sum  of  residues  at  simple  poles,  but  rather  includes  contributions  from 
branch  points,  the  branch  line  integral  is  ignored.  In  all  cases,  the  field  at  the  observation 
point  can  hence  be  written  as 

P{Robs\Rsp)  =  G{Robs\Rsp)~\' 

T  /  >  G j  (f)j  (^Zobs)[^j  {Robs  ■}  Rsp)  B j{Robs’)  -^sp)]  5 

^  i=l 

where 

Bj{Robs:  Rsp^  —  [  Bq  (^obs  —  shyJ^O  4^j{z^sh)^R{^sh')^sh\'^sp'j^sp^  ’  5 
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and  where 


5p)  /  P{Rsh\Rsp)y/kl  -  (  '^obs—shyJ^Q  ^j{^sh)^rdS  “h 

P{Rsh\Rsp)Ho''  (robs-sh\/ko  -  Aj^  4>'j{z,f,)aJS  . 

The  expression  Vobs-sh  is  the  radial  component  of  Robs  —  Rshi  and  the  a's  are  given  in  terms 
of  the  horizontal  radial  unit  vector,  the  vertical  depth  unit  vector,  and  the  target’s  outward 
surface  normal  unit  vector  as 


and 


Q^7>  —  Cf  *  Tlg}i 


-  ^Z  ' 


The  interpretation  given  to  the  above  equation  is  that  the  pressure  field  at  the  point 
Robs  consists  of  the  arrival  field  G,  which  has  not  been  scattered  by  the  target,  and  of  an 
arrival  scattered  by  the  target,  and  is  represented  by  the  linear  combinations  of  the  the  terms 
Bj  and  Aj. 
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APPARENT  TARGET  STRENGTH 


Knowledge  of  the  scattered  field,  the  source  strength,  the  transmission  loss  from  the 
source  to  the  phase  center  of  the  target,  and  the  transmission  loss  from  the  phase  center  of  the 
target  to  the  receiver  location  is  sufficient  to  calculate  an  “apparent  target  strength”  using 
the  sonar  equation.  The  qualifier  “apparent”  is  used  to  indicate  that  this  calculated  quantity 
will  be  different  from  the  comparable  quantity  in  a  free-field  environment.  This  quantity  will, 
in  fact,  depend  upon  the  specific  locations  of  the  source,  the  target,  and  the  observer.  This 
inability  to  treat  the  target  and  the  environment  separately  is  taken  to  be  a  violation  of  a 
fundamental  premise  of  the  sonar  equation.  The  formula  used  in  the  calculations  is 


T  S apparent  =  1 0  log 


( ^scattered  ®1  meter  Jrom  target 


V  h 


=  10  log 


incident  ©target 
^scattered  ©receiver 


incident  ®1  meter  from  source^ 


{TLsrc  to  TGT  +  TLtgj  to  RCVr)  , 


where 


and  where 


TLtgt 


to  RCVR 


10  log 


^  scatter  ed©receiver 


\ 


TL 


SRC  to  TGT 


10  log 


^ scatter ed©lmeter fromtar get  J 


incident  ©tar  get 


\  I incident©lmeter  from  source  ^ 
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INTERMODAL  SCATTERING  STRENGTH 


A  useful  way  to  look  at  the  problem  of  the  interaction  of  a  target  with  its  environment 
is  as  a  transition  operator  between  the  complex  amplitudes  of  the  ensonifying  and  scattered 
modes  of  the  duct.  This  is  a  more  detailed  description  of  the  scattering  process  than  is  the 
“apparent  target  strength,”  and  it  allows  for  a  deeper  understanding  of  the  phenomenon. 
In  particular,  this  formulation  of  the  problem  reveals  a  hidden,  approximate  relationship 
between  the  free-held  bistatic  target  strength  and  the  intermodal  transition  operator. 

Equation  (1)  can  be  rewritten  symbolically  as 

P{Rsg\Rsp)  =  Csrf  O  G{Rsp\Rsg)  ■ 

Then,  the  surface  pressure  can  be  rewritten  in  the  form 


-i(ro  tncVArg-Aj-7r/4)  /  , 


^  7rro  incy^O 


4^j{^sp)Dj(^Vsrf)  , 


where 


Dj{f,rf)  =  C,rf\e  Cine  ^>^0  j)\  , 


and  where  the  far-field  (Fraunhofer  zone  beyond  the  Fresnel  zone)  approximation 


'^srf—sp  —  ^0  inc  “f“ 


^srf  ’  ^0  inc 


has  been  applied. 

Similar  approximations  made  in  the  case  of  the  scattered  held  yield  the  following 
representation: 


P{Robs\Rsp)  —  G{Robs\Rsp)P 


EE 


o  o6s 'Y  ^0  Aj 


0  ^(^0 


j=i  fci  M  nro  obsJko  -  Xj  M  irro  inc-Jkl  -  A/ 


G jCi(l>j(^ZQ))g^<f>i(^Zgp')(f>j(^Zp^agg 

center  ase  center  )S,i  , 


where 


Sj,  =  [EA^o  inc')  ^0  065)  inc’)  se  center  ^4^l^^phase  center  j 


Eji{eo  ino  ^0  06s) 


^5r/‘^0  obs  /j  2 


Di{rsrf)e  ^  065 1 


^j4^j{_^sh)^r  “t"  4^j(^^sh}^z  dS  , 
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and  where 


me?  ^0  o6s)  — 


L 


^0  ob. 


-VW^) 


(l>j{Zsh)'^ Dt{fsrf)  ■  fishdS 


The  quantity  Sij  is  the  (linear,  rather  than  logarithmic)  analogue  of  the  free-held 
target  strength,  and  is  called  here  the  “intermodal  scattering  strength.”  It  is  the  extent  to 
which  this  quantity  is  not  a  constant  function  of  its  indices  that  determines  the  seriousness 
of  the  violations  of  the  sonar  equation,  since  the  sonar  equation  would  require  that  these 
terms  factor  completely  out  of  the  double  sum.  The  intermodal  scattering  strength  will  also 
be  a  symmetric  function  of  its  indices  by  virtue  of  acoustic  reciprocity  (the  symmetry  of  the 
Green’s  function).  Finally,  it  should  be  noted  that  the  intermodal  scattering  strength  is  a 
function  of  the  target  -physics^  the  target  depth  (generally)  and  the  propagation  modes  of  the 
environment.  While  the  intermodal  scattering  strength  matrix  need  not  neglect  the  Fresnel 
terms  (i.e.,  use  the  approximations  in  equation  (6))  in  the  ranges  between  the  target  and 
its  source  and  observer  points,  the  resulting  factorization  yields  a  form  for  equation  (7)  that 
is  computationally  efficient.  Specifically,  the  expensive  part  of  the  computations  (i.e.,  Sij) 
need  only  be  done  once  for  each  incident  target  aspect  and  observation  angle.  The  double 
sum  then  can  be  inexpensively  computed  for  any  combination  of  source  and  observer  range 
and  depth. 
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CALCULATED  RESULTS 


The  case  studies  of  a  rigid  sphere  and  an  elongated  rigid  body  in  simple  ducts  were 
chosen  to  demonstrate  the  model.  While  elastic  shells  are  easily  treated  with  this  methodol¬ 
ogy,  elasticity  adds  nothing  to  the  understanding  of  the  issues  being  addressed  in  the  present 
report.  Two  environments  were  investigated.  The  elongated  rigid  body  consisted  of  a  right- 
circular  cylinder,  with  a  length-to-diameter  ratio  of  6.3,  which  is  terminated  on  each  end 
by  a  hemisphere.  The  overall  length-to-diameter  ratio  of  the  body  is  thus  7.3.  The  first 
environment  to  be  investigated  was  a  channel  with  a  uniform  sound  speed,  a  uniform  depth, 
and  a  hard  bottom.  The  second  environment  was  a  channel  with  a  deep  sound  speed  axis 
that  gave  rise  to  convergence  zones.  The  sound  speed  profile  chosen  was  that  corresponding 
to  the  Monk  canonical  sound  channel  (Monk,1974). 

The  first  target/environment  configuration  to  be  considered  was  the  rigid  sphere  in 
the  uniform  duct.  The  center  of  the  sphere  was  on  the  depth  centerline  of  the  duct  at  a 
distance  of  1000a  from  the  source  and  ka  «  1.  First,  the  scattered  field  is  calculated  along 
a  locus  of  depths  from  the  surface  to  the  channel  bottom  and  lying  at  the  same  range  as 
the  source  (i.e.,  along  a  vertical  line  containing  the  source  point  as  detailed  in  figure  1).  It 
is  instructive  to  compare  the  apparent  target  strength  calculated  from  the  scattered  field 
from  the  rigid  sphere  in  an  unbounded,  uniform  environment  with  that  for  the  sphere  in 
a  bounded,  uniform  channel  (figure  2)  and  also  with  that  in  the  convergence  zone  (figure 
3).  Clearly,  the  effect  of  the  environmental  coupling  is  stronger  for  the  case  of  the  shallow 
environment  with  a  hard  bottom  than  for  the  convergence  zone.  This  effect  is  believed  to  be 
a  consequence  of  the  different  ranges  of  modal  phase  velocities  (and,  hence,  effective  incident 
angles)  for  the  two  environments.  Specifically,  in  the  case  of  the  convergence  zone,  the  phase 
velocities  associated  with  the  modes  carrying  most  of  the  energy  into  and  away  from  the 
target  correspond  to  angles  within  ±1  degree  of  horizontal.  In  the  case  of  the  shallow  duct, 
on  the  other  hand,  the  modal  phase  velocities  correspond  to  much  greater  angles.  It  is  this 
essentially  bistatic  scattering  phenomenon  that  gives  rise  to  the  large  differences  between 
the  apparent  target  strength  in  the  shallow  duct  and  the  free-field  target  strength. 

As  a  further  example,  the  apparent  target  strength  of  the  elongated  body  in  the  shallow 
duct  is  shown  in  figure  4  for  the  case  kaK,  1.  The  target  sits  at  end  incidence  to  the  direction 
of  ensonfication  and  observation.  It  is  noteworthy  that  the  apparent  target  strength  of  the 
elongated  object  shows  even  stronger  excursions  from  the  predictions  of  the  sonar  equation 
than  does  the  sphere. 

Examples  of  the  importance  of  multiple  scattering  is  shown  in  figures  5  through  8. 
In  figures  5  and  6,  the  apparent  target  strength  is  shown  for  the  case  of  the  sphere  located 
at  distances  of  7a  and  2a  from  the  surface,  respectively.  The  effect  of  multiple  scattering 
was  determined  to  be  negligible  at  the  center,  but  more  significant  for  the  target  near  the 
surface.  Corresponding  apparent  target  strengths  are  shown  in  figures  7  and  8  for  the  target 
near  the  rigid  bottom. 
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The  effects  of  neglecting  the  Fresnel  terms  implicit  in  the  use  of  the  approximation 
of  equation  (6)  are  shown  in  figure  9.  The  case  of  the  elongated  body  oriented  at  beam 
aspect  to  the  source/receiver  axis  is  treated.  A  comparison  of  the  complete  solution  with  the 
approximate  solution  indicates  that  the  effect  of  being  in  the  Fresnel  zone  is  a  second-order 
effect,  but  is  clearly  detectable  at  moderate  ranges. 
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CONCLUSIONS 


Numerical  confirmation  has  been  obtained  of  the  significant  violations  of  the  sonar 
equation  which  can  occur  in  dispersive  environments.  This  problem  is  shown  to  be  more 
serious  for  a  shallow  water  duct  having  a  hard  bottom  than  for  a  convergence  zone  envi¬ 
ronment.  A  simple  interpretation  of  this  phenomenon  is  that  the  great  range  modal  phase 
velocities  in  a  shallow  duct  gives  rise  to  stronger  equivalent  bistatic  effects. 

Multiple  scattering  effects  have  been  shown  to  be  observable  for  a  target  near  the 
surface  of  the  duct.  They  were  not  significant  near  the  center  of  the  duct  under  study. 

The  effects  of  neglecting  Fresnel  terms,  which  makes  the  problem  particularly  efficient 
to  solve,  are  not  entirely  negligible  at  useful  ranges. 
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14.6  a 
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Figure  1 .  Source/target/receiver  arrangement  for  shallow  duct. 
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Figure  2.  Apparent  target  strength  of  hard  sphere  in  uniform  duct  with  hard  bottom. 


Figure  3.  Apparent  target  strength  of  hard  sphere  in  convergence  zone. 
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Figure  6.  Multiple  scattering  effects  on  apparent  target  strength  of  rigid  sphere  near 
surface.  Target  Depth/a  =  2. 
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Figure  7.  Multiple  scattering  effects  on  apparent  target  strength  for  rigid  sphere  near 
rigid  bottom.  Target  Depth/a  =  53. 
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Figure  8.  Multiple  scattering  effects  on  apparent  target  strength  for  rigid  sphere  near  hard 
bottom.  Target  Depth/a  =  58. 


Range/a 

Figure  9.  Effect  of  Fresnel  terms  on  apparent  target  strength  for  rigid  cylinder  with 
end  caps.  Target  at  center  of  duct  and  at  beam  aspect. 
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APPENDIX 

EQUATION  DERIVATIONS 


A-l 


A.l  FORMULATION  OF  THE  PROBLEM 


This  appendix  is  a  generalized  version  of  equation  derivations  detailed  by  McDaid  et 
al.  (1992).  The  general  approach  in  solving  this  problem  is  to  use  a  version  of  the  Helmholtz 
Integral  Equation  that  has  been  modified  to  account  for  the  refractive  acoustic  environment 
in  combination  with  a  so-called  normal  mode  formulation  of  the  acoustic  propagation.  The 
refractive  Green’s  function  is  used  to  calculate  the  pressure  incident  on  the  wetted  surface 
of  the  shell.  The  baseline,  free-held  version  of  the  Helmholtz  Integral  Equation  is  then  used 
to  calculate  the  pressures  and  velocities  on  the  wetted  surface  (equation  A. 6).  Finally,  the 
refractive  version  of  the  Helmholtz  Integral  Equation  is  used  to  calculate  the  far-held  scat¬ 
tered  pressure  (equation  A. 5).  The  modal  representation  of  the  refractive  Green’s  function 
(equation  A. 9)  makes  the  problem  particularly  easy  to  formulate. 

A.2  EXTERIOR  HELMHOLTZ  EQUATION  FOR  REFRACTIVE  ENVIRON¬ 
MENT 

The  formulation  of  the  scattering  problem  for  a  shell  in  a  refractive  environment  in 
which  the  sound  speed  prohle  is  horizontally  stratihed  is  straightforward,  since  it  is  so  similar 
to  the  equivalent  problem  in  a  uniform  acoustic  medium.  The  derivation  of  the  relevant 
equations  is  included  herein  for  the  sake  of  completeness,  rather  than  as  a  demonstration  of 
their  novelty.  Throughout  this  work,  the  time  dependence  of  the  signals  is  assumed  to  be 
The  domain  of  interest  is  a  half  space,  V,  with  a  pressure  release  boundary  condition 
at  the  surface.  Supper-  Let  the  vectors  Rfg  and  Rgg  be  the  field  and  source  points  for  the 
Green’s  function  for  the  following  boundary  value  problem 

G{Rjg\Rsg)  =  -S{Rfg  -  Rsg)  in  V  , 

G{Rjg\R  sg)  —  0  on  Supper  ? 

and 

T  I  -7^1  n 

lim  r.  +  ikG\  ~  0, 

Tg-^oo  dr 

where  Vg  =  \Rjg  —  and  k[zfg)  ~  ujlc{zjg).  The  term  Ufi  is  the  outward  normal  on  the 
fluid  body  at  the  shell  surface,  Ssh-  Note  that  nji  =  —fish-)  where  Ush  is  the  outward  normal 
to  the  elastic  shell,  and 

g~^sg\ 

lim^  G{Rfg\R,g)  -  — ^ - g—  =  0. 

Rfg-*  Rsg  [  4:'K\Rfg  -  Rsg\ 

This  last  equation  is  used  to  justify  the  use  of  the  free-space  Green’s  function  in  equation 
(A.6)  as  an  approximation  of  G.  Also  note  that  the  Green’s  function  is  symmetric,  i.e., 
G{Rjg\R,g)  =  G{Rsg\Rig). 

Let  the  vectors  Rjp  and  Rgp  be  the  field  and  source  points  for  the  pressure  in  the 
following  boundary  value  problem 

[^Rfg  +  k'^i^fp)]  P{Rfp\Rsp)  =  0  in  U  , 
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P{Rfp\R 

sp)  —  0  on  Supper  , 


^ RjpP{Rjp\Rsp)  ■  nji  =  T[P{Rfp\Rsp)]  on  Ssh  , 


and 

,dP 

lim  rp\— — h  ikP\  =  0  , 

Tp^oo  (jr 

where  Vp  =  \Rfp  —  i^spl,  k{zfp)  =  ulc{zfp)^  and  Ssh  is  the  wetted  surface  of  the  elastic  shell. 
The  functional  relationship  P’  is  used  to  embody  the  effect  of  the  elastic  shell.  Note  that  P 
is  a  mapping  between  functions. 

In  the  present  example,  let  Rfp  —  Rfg  =  Rj.  We  then  have  the  result 


P(.R,\R.r)  [V|,  +  G(R,\R„)  -  G(R,\R.,)  [vjj,  +  e(z,)]  P{R,\R,,)  =  0  , 


or 

P(R,\R.,)Vl^G{R,\R.,)  -  G(R,\R„)V%P{Rj\R„)  =  0  .  (A.l) 

The  following  identities  are  useful  in  simplifying  the  foregoing  expression: 

V„,  .  [P{Rf\R,,)Vs,G[R,\R.,)]  =  Vk,P{R,\R.p)  ■  r,G{R,\R.,)  + 

P{R,\R.,)V},,G(R,\R.,)  =  0  , 

and 


Vr,  ■  [G(RAR.,)'^r,P(R,\R„)]  =  \!r,G{R,\R..)  ■  Vr,P{R,\R.p)  + 

G(Rj\R.,)V%P{RAR„)  =  0  . 

Hence,  equation(A.l)  can  be  written  in  a  simplified  form  as 

Vh,  •  [P{Rf\Rsp)VR,GiRf\Rs,)  -  G{Rf\Rsp)VR^PiRj\Rsp)]  =  0. 


This  expression  can  be  integrated  over  the  volume,  which  excludes  the  shell  and  its 
interior,  a  tiny  sphere  of  radius  e  centered  at  R^p,  and  another  tiny  sphere  or  radius  e  centered 
at  Rsg.  Since  there  are  no  sources  in  this  volume,  one  has  the  result 

Lr  „  ^R,-lP(Rl\R.p)^n,G(Rl\R»)- 

—  Keg  — Kep— Kg/i 

G{Rf\Rsg)VR^PiRf\Rsp)]dV  =  0. 

This  is  readily  converted  into  a  surface  integral  of  the  form 

/  [P{Rf\Rsp)VR,GiRf\Rsg)-G{Rf\Rsg)VR^P{Rj\Rsp)]-nfidS  =  0.  (A.2) 

J  S^g+Sep+Ssh 
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The  integrals  over  the  surfaces  S^g  and  Sf^p  have  particularly  simple  limiting  forms.  Note 
that 


/  P{Rf\R,p)VngG{RjAg)-njidS 

Seg 


PiRjR^p)  [  VR^GiRf\R,,)-njidS 

Scg 


P{Rsg\Rsp) 

P{Rsg\Rsp) 

P{Rsg\Rsp) 


^-ikc 

47re 


) 


—  (47re^) 


( 

\  47re 

47re^  j 

Similarly,  one  has 

/  G{Rs\R,g)Vn,P{Rj\R,p)-njidS 

J  S^p 


Note  also  that 

/  G{Rf\R,g)VR^P{Rj\R,p)  ■  npdS 

JScg 


Similarly,  one  has 

/  P{Rf\R,p)VR,G{Rf\Rsg)  ■  npdS 

Sep 


G{R,p\Rg,)  (  Wr^P{Rj\R,p)  ■  njidS 

Sen 


~ikc 


G{Rsp\Rsg) 
G{Rsp\Rsg)  ■ 


—ike  ^—ike\ 


-(A'kc)  —ik - - 

\  47re  47re^  J 


\VR^P{Rsg\Rsp)\  [  G{Rf\Rsg)e-nfidS 

Sep 

\Vr^P{R,p\R,p)\^  f  e-npdS 

4r  TT  C  1/  Sep 

0. 

[  P{Rf\R,p)e-nf,dS 

s eg 

|VR^G(.R5p|i?5p)|- —  /  e-UfidS 
i  TT  C  «/  Seg 

0. 


These  four  approximations  are  exact  in  the  limit  as  e  shrinks  to  0.  The  surface  integral  can 
thus  be  written  as 


P{Rsg\Rsp)  ~  G[Rsp\Rsg)p 

f  [P{RsH\RspWR,,G{Rsh\Rsg)  -  G{Rsh\Rsg)VR^,P{Rsk\Rsp)]  •  njidS  =  0.  (A.3) 

■’Ssh 
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The  surface  integral  can  thus  be  rewritten  in  terms  of  the  shell  normal  as 
P{Rsg\Rsp)  =  G{Rsp\Rsg)  + 

[  [PiRsh\Rsr)^R,MRsh\Rsg)  -  GiRskA9)^R.,P{Rsh\Rsp)]  •  nshdS 
Js.h 


(A.4) 


A  change  of  notation  can  be  invoked  and  the  symmetry  of  the  Green’s  function  can 
be  used  to  cast  the  problem  into  a  form  familiar  to  those  working  in  scattering  theory,  such 
that 

P{Robs\Rsp)  =  G{Robs\Rsp)+ 

'^Ssh 

where  R^bs  is  substituted  for  the  term  Rgg. 

The  interpretation  of  the  above  result  is  that  the  resulting  field  consists  of  the  direct 
arrival  from  the  source,  represented  by  G,  and  a  term  scattered  off  the  shell  and  represented 
by  the  surface  integrals  over  the  shell  surface. 

A.3  EXTERIOR  HELMHOLTZ  EQUATION  WITH  REFRACTIVE  ENVIRON¬ 
MENT  AND  LIQUID  BOTTOM 

In  the  case  of  a  “liquid  bottom,”  one  has  essentially  two  different  media  with  which 
to  deal.  It  is  most  correct  to  specify  the  propagation  in  terms  of  two  Green’s  functions,  one 
for  each  of  the  media.  The  Pekeris  duct  is  a  special  case  of  this  situation.  In  this  section,  we 
will  use  the  convention  of  calling  these  functions  Gi  and  G2  for  the  Green’s  functions  in  the 
surface  layer  and  basement  layer,  respectively.  The  analysis  for  this  case  follows  that  of  the 
single  medium  very  closely.  The  difference  is  that  boundary  conditions  at  the  liquid-liquid 
interface  must  be  satisfied.  These  results  are  easily  extended  to  the  case  of  a  multilayer 
liquid  layer. 

The  target  will  always  be  assumed  to  lie  in  the  surface  layer.  An  equation  correspond¬ 
ing  to  equation  (A. 2),  but  applying  in  the  surface  layer,  is 


[P,{Rf\R,,)VR^G^iRj\R,g)- 


Gr{Rf\Rsg)WR^Pr{Rf\Rsp)]  ■  dfndS  =  0. 

The  subscript  fll  refers  to  the  outward  normal  on  the  surface  layer  at  the  interface,  and  the 
subscript  b  refers  to  the  interface.  Similarly,  one  has,  for  the  field  in  the  basement  layer,  the 
equation 


[P2{Rf\Rsp)VR,G2{Rf\Rsg)- 


G2{Rf\Rsg)^R,P2{Rf\Rsp)]  ■  df,2dS  =  0. 
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The  subscript  fl2  refers  to  the  outward  normal  on  the  basement  layer.  In  the  foregoing 
equations,  provision  has  been  made  for  the  possibility  of  the  source  and  field  points  lying  in 
either  of  the  two  layers. 

The  two  media  are  coupled  through  the  acoustic  boundary  conditions  which  obtain  at 
the  interface,  i.e.,  the  equality  of  stress 

G,{Rj\R,,)  =  G,{RfA,)  , 
and  the  equality  of  normal  velocities 

RfGi{Rf\Rsg)  ■  nfti/ Pi  -  -V ji^G2iRf\Rsg)  •  nji2/ p2  , 
where  it  is  noted  that  njn  =  — n//2. 

Similarly,  the  pressure  fields  satisfy  the  same  boundary  conditions 

Pl{Rj\R,g)  =  P2{RsAg) 

and  the  equality  of  normal  velocities 

^ RfP\{R}\Rsg)  •  Rflllpi  =  —VjijP2{Rf\Rsg)  ■  nfi2l P2  • 


A  complete  treatment  of  the  sound  field  would  require  that  a  number  of  cases  be 
considered,  in  accordance  with  the  layers  in  which  the  source  and  field  points  lie.  In  this 
study,  the  most  important  case  (and  the  only  one  for  which  calculations  will  be  done)  is  that 
one  in  which  both  Rgg  and  R^p  lie  in  the  surface  layer. 

The  equations  for  the  pressures  in  the  two  media  become  as  follows: 


Pl{Rsg\Rsp)  ~  f^l(-Rsp|.Rsg)  + 

/  [Pi{Rsk\Rsp)yR^,,Gi{R,f,\R,g)  -  Gi{R,i,\R,g)Vn^^Pi{R,t,\R,p)]  -  + 

j  [Pi{Ri\Rsp)VR,Gi{R,\R,g)  -  Gi{R,\R,g)Vn^Pi{R,\R,p)\  •  n^ndS  =  0  , 

J  big 

j^[P2{Rb\RspWRfi2{Rb\Rsg)  -  G2{Rb\Rsg)VR,P2{Rb\Rsp)]  '  njl2dS  =  0. 
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Using  the  interface  boundary  conditions,  one  has 


0=  / 

J  Si 

—  [  [Pl{^b\Rsp)'^ Rt, 

JSi 


GiiRblRsg)  -  Gi{Ri,\Rsg) 


~P2 

Pi 


^ Rt,RiiRb\Rsp)]  ■  nfi2dS 


-P2 


pi  JSt 


f  [PiiRb\Rsp)^R,G^{R,\R,,)  -  G^{Rb\Rss)VR,P^{RbA,)]  ■  njndS 

JSh 


If  one  substitutes  this  latter  result  into  the  equation  for  Pi ,  then  one  has  the  expected 
result,  which  is  equivalent  to  equation(A.3). 

PliRsg\Rsp)  ~  Gi{Rsp\Rsg)-\- 

f  [Pl{Rsh\Rsp)yR,,Gi{R,h\Rsg)  -  Gi{R,h\Rsg)^R,,Pl{Rsh\Rsp)]  •  njidS  =  0. 

'^Sgh 


By  using  the  appropriate  change  of  notation  and  invoking  the  symmetry  of  the  Green’s 
function,  one  arrives  at  an  equation  identical  to  equation(A.5). 

A.4  SURFACE  INTEGRAL  EQUATION 

In  this  case,  perform  the  same  analysis,  except  that  one  lets  Rgg  be  on  the  surface  of 
the  shell.  In  this  case,  one  has  the  equation 


-P(Rsg\Rsp)  —  G{Rsp\Rsg)-\- 


[  [PiRshAp)^R.HG{Rsh\Rsg)  -  G{Rsk\Rsg)^R,,P{Rsh\Rsp)]  ■  n.ndS  . 

'^Ssh 


It  is  this  equation  that  must  be  solved  in  order  to  evaluate  the  pressure  and  velocity  on 
the  surface  of  the  shell.  In  the  present  instance,  the  approximation  is  used  in  which  the  free- 
field  propagation  Green’s  function  is  used  within  the  integral,  rather  than  the  actual  refractive 
Green’s  function.  The  forcing  function,  G{Rsp\Rsg),  for  the  integral  equation  is,  however, 
the  actual  refractive  Green’s  function.  More  will  be  said  regarding  this  approximation  in 
section  A.ll.  For  the  sake  of  clarity,  this  integral  equation  can  be  rewritten  in  terms  of  the 
free- field  Green’s  function  Go  as 

-P{Rsg\Rsp)  =  G{Rsp\Rsg)-{- 

I  [P{Rsh\Rsp)'^ R,hGo{Rsh\Rsg)  —  Go{Rsh\Rsg)'^ R.^PiRshlRsp)]  ‘  flshdS  .  (A. 6) 

'^Ssh 
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If  one  denotes  the  inverse  of  the  foregoing  integral  operator  by  the  symbol  £srf,  then 


p(ii..Ap)=  I  K{Ap,Ap;R,,)G{R,Mp)<is(A'j, 

or 

P{Rsa\Rsp)  -  C,rjG{Rs^\R,g)  . 

A.5  NORMAL  MODE  REPRESENTATION  OF  PROPAGATION  GREEN’S 
FUNCTION 

In  the  present  study,  two  categories  of  stratified  sound  speed  profiles  have  been  con¬ 
sidered.  In  the  first  instance,  an  approximate  solution  has  been  found  for  the  so-called 
“Pekeris”  wave  guide.  This  environment  consists  of  two  fluid  layers,  each  with  a  uniform 
density  and  sound  speed.  The  upper  layer  is  finite  in  depth  and  is  bounded  above  by  a 
pressure  release  boundary.  The  lower  layer  is  an  infinite  half  space.  In  the  second  instance,  a 
multilayer  environment  is  chosen  in  which  propagation  exhibits  convergence  zone  behavior. 
The  sound  speed  is  a  continuous  function  of  depth.  Within  each  layer,  the  square  of  the 
sound  speed  has  a  hyperbolic  dependence  on  depth.  In  the  deepest  layer,  which  is  an  infinite 
half-space,  the  sound  speed  has  an  asymptotic  limiting  value  of  zero. 

The  propagation  Green’s  function  is  a  solution  to  the  following  boundary  value  problem 

[vi„  +  G(Ri,A,)  =  -S{Rj,  -  R„)  in  V  , 

G(RfgjRsg)  =  0  on  Supper  , 
dG 

where  Vg  is  the  horizontal  separation  distance  between  Rfg  and  Rsg.  In  the  case  of  a  horizon¬ 
tally  stratified  sound  speed  profile,  it  is  convenient  to  solve  the  sound  propagation  problem 
in  cylindrical  coordinates.  In  those  cases  where  the  Green’s  function  can  be  represented  as 
a  series  of  eigenfunction  products  or  a  residue  series,  without  the  inclusion  of  a  branch  line 
integral  contribution,  one  has 

'  CO 

G{rj,zf\r,,z,)  =  Xg)G,{zf\z,-,-\g)  , 

^  j=i 

where  the  summation  is  over  a  discrete  set  of  values  of  the  (complex)  separation  parameter 
A,  where  r  is  the  horizontal  separation  between  Rj  and  Rs,  and  where  ko  is  a  reference 
wave  number.  The  choice  of  a  value  for  ko  is  dependent  on  the  sound  speed  profile  being 
considered  in  a  particular  application. 


A-8 


In  the  present  work,  a  standard  notation  will  always  be  used  when  defining  cylindrical 
coordinate  systems.  The  depth  coordinate  ^  will  always  be  positive  downward,  so  will 
point  down.  The  radial  coordinate  r  will  be  horizontal,  and  the  following  subscripting 
notation  will  be  used  for  r  and  the  horizontal  radial  unit  vector  e^: 


rj-s  =  \{Rf  -  Rs)  -  [{Rf  -  Rs)  ■  (A. 7) 

and 

Or  =  \{Rf  -  Rs)  -  [{Rf  -  Rs)  ■  e^]e^|/r/_s  •  (A. 8) 


The  depth  Green’s  function,  is  further  written  as 

Gz{^f\^s]  ~^j)  “  Gj(f)j[zf)(l)j[Zs)  , 

where  the  function  (j>j{z)  is  the  j-th  depth  “eigenfunction.”  The  reason  for  temporizing 
with  regard  to  the  name  of  the  function  <f)j  is  that  the  series  is  actually  a  residue  series 
in  the  cases  of  certain  sound  speed  profiles,  and  it  is  not  clear  if  the  functions  should  be 
considered  an  orthonormal  basis  for  the  function  space  under  consideration  since  they  may 
not  be  integrable.  The  functions  (j)j  have  been  chosen  such  that  =  0  and  d(f)jldz  —  1 
for  2:  =  0,  where  the  depth  2:  is  a  positive  quantity  and  has  a  value  of  0  at  the  surface. 

By  substituting  this  result  into  the  original  equation  for  the  Green’s  function,  one  has 


•  00 

G{rf^  ZfVs,  Zs)  =  \^  -  \j)(l)j{zf)4>j{zs)  .  (A.9) 

^  i=i 

The  gradient  of  the  Green’s  function  is  given  as 

^  i  °° 

VG{rf,Zf\rs,Zs)  =  jY.Gj 
^  i=i 

[er\/kl  -  \jHP  [ryjkl  -  4)f{zf)(l)j{zs)  + 

.  (A.IO) 

The  horizontal  radial  unit  vector  Cr  points  away  from  the  source  point,  and  the  axial 
unit  vector  points  down,  since  the  sign  convention  will  be  that  2:  is  depth,  and  hence, 
increases  in  the  downward  direction  (equations  A. 7  and  A. 8). 

The  expression  for  the  surface  pressure  can  hence  be  rewritten  as 

P{Rsrf\Rsp)  =  4  Z!  GjCsrf  \^Ho\rsrf-3p\/kl  -  Xj)<f>j{Zsrf)j  M^sp)  • 
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Using  the  asymptotic  representation  of  the  Hankel  function,  one  has 


4'ji^sp) 


Let  the  location  of  the  phase  center  of  the  scatterer  be  denoted  Rphase  center  and  let 


ro 


=  R 


phase  center 


R 


sp 


se  center 


be  the  vector  from  the  source  point  to  the  phase  center  of  the  target.  Let 


r  f 

srj 


=  Rgrf  —  R 


phase  center 


'srf 


R 


■phase  center 


)  ■ 


be  the  vector  from  the  phase  center  of  the  scatterer  to  the  surface  point.  The  following 
approximation  is  used  for  large  values  of  the  horizontal  range  from  the  source  point  to  the 
surface  point,  r^rj-spi 


^srf—sp  —  ^0  inc  T 


^ srj  '  ^0  inc 

j  •> 

r  0  inc 


where  Tq  inc  is  the  vector  pointing  from  the  source  point  to  the  phase  center  of  the  target, 
and  ro  inc  =  inc\-  The  pressure  can  be  written  approximately  as 


P{Rsrs\Rsp)  = 


A  ^  CjCsrf 
^  i=i 


IN 


^^0  ir 


'\J^0  ~ 


-iiro  i  <f>jiZsp)  • 


An  approximate  “modal  factorization”  can  be  performed  as  follows: 


P{Rsrf\Rsp) 

i  "" 


_ 

i=l  ^  TTCq  mcy  ^0 


■'srf 


■7^V/'^0  >nc 
e  ^  ^0  inc 


^^UjiZsrj) 
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or 


P(R.rAR.p)  =  j'£Cj 

^  i=i 


N 


ir 


■“’■•  ,  (A.n) 


where 


Di(Ptrf)  =  C.„ile  ’•  'C,.,  \/*"  ^>Vj(2.r/) 


(A.12) 


Equation  (A. 11)  is  simply  the  surface  pressure  that  has  been  calculated  by  using  CHIEF /MARTSAM 
for  the  case  of  a  point  source  located  at  Rsp.  Equation  (A.12),  on  the  other  hand,  is  the 
surface  pressure  that  would  be  calculated  by  using  the  partial  point  source  strength  that  is 
given  by 

e  ’•0  .nc  V  ^’(j,-[Zsrf)  , 

rather  than  the  entire  source  strength  G{Rsrf\Rs-p)- 

Similar  results  can  be  obtained  for  the  gradient  of  the  surface  pressure, i.e., 


VP{R^rj\Rep)  = 


4  ^ 

^  j=i 


or 


^P(R.„\R.p)  =  2Y:c, 

^  j=l  ^ 


Trro  ifi 


The  errors  in  this  expression  are  C>(l/ro  inc)- 

A.6  MODAL  REPRESENTATION  OF  SCATTERED  PRESSURE 

By  substituting  equation  (A. 9)  and  equation  (A. 10)  in  equation  (A. 5),  the  two  surface 
integrals  can  now  be  evaluated  in  terms  of  quadratures  involving  the  depth  eigenfunctions. 
In  particular,  one  has 

J  G{Robs\Rsh)'^ R,hPi^sh\Rsp)  ■  n^hdS  =  ^  H G j(t>j{zohs)Bj{Robs-, Rsp) , 

Ssh  4 


where 


R j  i^Robs )  R  5p)  ~sh  \J~^0  R {j' sh')  ^ sh\^ sp')  ^ sp^  ’  , 

•^Ssh 
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and  Tobs-sh  is  defined  in  accordance  with  the  conventions  of  equations  (A. 7)  and  (A. 8). 


Similarly,  one  has 


/  P{Rsh\Rsp)'^ R,HG{Robs\Rsh)  ■  UshdS  =  -  J2^J^jiZobs)Aj{Robs,Rsp)  , 


where 


Aj{Robs,  Rsp)  =  [  P{Rsh\Rsp)\/ko  -  \jH'^^\robs-shJ  kl  -  \j)(j)j{zsh)oLrdS+ 

I  P{Rsh\Rsp)Ho\robs-shA/ko  -  Xj)(l)'j{zsh)cX;,dS  , 

and  where  the  a's  are  given  in  terms  of  the  unit  horizontal  radial  vector  and  the  unit  vertical 
depth  vector  as 


and 


CXr  6y.  *  Tl , 


• 


It  is  clear  from  equations  (A. 4)  and  (A. 10)  that  the  radial  unit  vector  must  point  away 
from  Robs- 

The  following  asymptotic  approximations  of  the  Hankel  function  and  its  derivative  are 
used  in  the  calculations  (equations  9.2.4  and  9.1.28  of  Abramowitz  &  Stegun  (1964)) 


TTX 


and 


~  -i\[^ 

V  TTX 


The  field  at  the  observation  point  can,  hence,  be  written  as 


PiRo6s\Rsp)  =  GiRoBs\Rsp)  + 

i  ^ 

T  ^j^i('^o65) [^j (-Ro6s?  Rsp)  Bji^Robs’i  ■  (A. 13) 

i=l 

The  interpretation  given  to  the  above  equation  is  that  the  pressure  field  at  the  point  Robs 
consists  of  the  arrival  field  G,  which  has  not  been  scattered  by  the  target,  and  of  an  arrival 
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scattered  by  the  shell  and  represented  by  the  linear  combinations  of  the  the  terms  Bj  and 

A,. 


Using  the  approximate  modal  factorization,  one  has 


be  the  vector  from  the  field  point,  where  the  scattered  pressure  is  being  calculated,  to  the 
phase  center  of  the  scatterer  .  Use  the  approximation 


where  Tq  obs  =  i^h  o6s|)  to  get  the  result 
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This  result  can  be  written  with  simplified  notation  as 


where 


inc^  ^0  obs}  — 


Similar  approximations  can  be  found  for  the  terms  B, 
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or 


where 

Fjl{€.Q  iYlO  ^0  o6s) 

/  •  H,hdS  . 

Js,h 

The  scattered  pressure  (equation  A.  13)  can  be  rewritten  as 
P{Robs\Rsp)  =  G{Robs\Rsp)+ 


Cj  Cl<j)j{Zobs)4’l{  ^sp  )  [^jl  (  ^0  mc)  ^0  o6s)  Fjl{^0  ino  ^0  o6s)]  • 


It  is  the  quantities  Ejiie^  inc-,^0  ohs)  Fji{eQ  inc^^o  ohs)  that  are  the  fundamental 
scattering  coefficients  which  must  be  calculated.  The  analogue  of  a  target  strength,  call  it 
the  “cross-modal  target  strength,”  or  the  “inter-modal  target  strength,”  can  be  written  as 

TSjl  —  [£'j/(eo  me?  ^0  obs)  ~  Fjl{eQ  mc?^0  obs)]/[^j{^ phase  center) 4^l{^phase  center)]  • 

The  equation  for  the  scattered  pressure  can  thus  be  written  as 
P{Robs\Rsp)  =  G{Robs\K)  + 


C jCl(l)j(^ZQl)s)(j>l{Zsp)(j)j(^Zptiase  center) 4^l{^phase  center)T S jl 
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A.7  NUMERICAL  APPROXIMATIONS 

In  the  version  of  the  CHIEF/MARTSAM  software  to  be  used,  an  approximation  is  to 
be  used  in  which  surface  elements  are  assumed  to  have  constant  pressure  and  normal 
velocity.  In  this  case,  the  coefficients,  Bj  and  Aj,  appearing  in  equation  (A.  13)  have  the 
following  representations: 


B jn{_Robs'i  Ftsp)  —  / 

Js 


( 7* 065-5/1  \/^-Aj)-7r/4 


TT 


Pij' shi  sp"}  ^sp)  '  • 

This  can  be  rewritten  approximately  as 


P jni^Pohs')  Psp  )  Pij'shn’)  '^sp’i^sp)  ’  '^shn  ^  ^ 


L 


■n  ^ 


(^^obs—shyJ^O  ) 


-i^PnVnB'^/'^Qyjn  ■ 


The  total  expression  for  Bj  is 


B,{Robs,  Rsp)  =  J2  PnVn{Rsp)Q,jn{Robs) 


N 


n=l 


In  the  case  of  the  coefficients  Aj,  one  has 
^jn  (  Robs  ?  Rsp  )  ~ 

Is,  P{R.k\R.,)^K  -  A,  (- 


(To65  -  56  \AfT^)  -  TT /4 


fsn  R(Rsh\Rsp) 


This  can  be  rewritten  as 


'^(robs-sh\/k^-^j) 


(»’o6«-s/i\Af^A7)-7r/4j 


>jizsh)ardS+ 


<i>'j{Zsh)0CzdS  . 


AjniRobs^Rsp)  =  I 

JSn 

[-i^kl  -  Xj(t)j{Zsh)oir  +  <i>'j{Zsh)ci^]dS  . 


{p'obs—sh'\J  Rq  Aj) 
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As  before,  this  again  can  be  rewritten  in  the  approximate  form 


A,n{Robs,Rsp)  =  P{Rsh\Rsp)e^^^^  I 

J  Sn 


ohs~sh\J^Q 


=  PnP^'^Q 


dS 


'ypjn  • 


The  total  expression  for  Aj  is 


Aj{Robs,Rsp)  =  P"^^ 


N 

Y.Pn(R  sp^Q  pjni^Robs)  * 

71=1 


There  is  a  rough  correspondence  between  the  coefficients  Bj{Robs^  Rsp)  Aj{Robs^  Rsp) 
and  the  coefficients  Bjj  and  Ay/  as  used  in  the  formulation  of  the  baseline  CHIEF/M  ARTS  AM 
software.  That  correspondence  is  as  follows: 


B jiyRohs'i  Rsp)  "} 

and 

(  Robs  ?  R  sp 

A  noteworthy  difference  between  the  baseline  version  of  the  CHIEF /MARTS AM  soft¬ 
ware  and  the  case  of  the  refractive  environment  is  that  in  the  latter  case,  one  cannot  extract 
a  simple  factor  which  will  account  for  the  propagation  loss  of  the  scattered  wave  away  from 
the  shell. 

The  expression  for  the  total  resulting  field  generated  by  the  source  in  the  presence  of 
boundaries  and  the  shell  is  hence  written  as 


P{Robs\Rsp)  =  G{Robs\R  sp)  P  P scat{Robs\Rsp) 


where  the  contribution  of  the  scattering  from  the  shell, is  written  as 


Pscat{Robs\Rsp)  =  JP^^"J2GjM^obs) 
^  i-1 


N  N 

sp  )Qpjn{Robs)  “I"  ^  V  Pn^n{Rsp)Qvjn{Robs) 


,n=l 


71  =  1 
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The  order  of  summation  (in  n  and  j)  has  been  written  arbitrarily  in  the  foregoing 
equation.  An  issue  that  must  be  resolved  in  the  course  of  the  study  is  the  determination  of 
which  order  of  summation  is  more  computationally  efficient. 

The  quadratures  can,  in  turn,  be  represented  as  weighted  sums  of  values  of  the  in¬ 
tegrands,  since  a  Gaussian  quadrature  procedure  will  be  used.  The  terms  can  hence  be 
rewritten  as 


Q  vjn 


Is 


'^{j^ohs—sh\J^^ 


(^o6s -.h  \/^  >^3  )  )  dS 


Mn 
j=l 


zjnm  •> 


where  is  a  set  of  Gaussian  quadrature  weights,  and  where 


rjnm 


and 


\  T^{rohs-sh\jH  -  Xj) 

fzjnm  —  • 


For  the  sake  of  simplicity,  the  obvious  dependence  of  robs-sh  and  Zsh  on  the  indices  m 
and  n  has  been  suppressed.  Similarly,  one  has  that 


Qpjn  =  [  e  ^(^obs~shy/kl  Xj)  _ 

\  ^{rohs-sh\lkl  -  A^) 

[  ^  \/^0  ( "^5/1 ) “1“  4^j{^sh)^z]dS 


Mn 

i=l 


rzjnm  ? 


where 


f rzjnm  ^^^ji^sh}^r  "b  ^ j^^sh)^z  • 


A-18 


A.8  PEKERIS  WAVEGUIDE 


In  the  case  of  this  environment,  a  surface  liquid  layer  of  depth  h  lies  over  a  liquid 
half-space  having  different  acoustic  properties.  The  specific  form  of  the  sound  speed  profile 
chosen  as  a  function  of  depth,  z,  is 


,if  0  <  z  <  h 
I ,if  /j  <  z 


The  density  is  also  nonuniform,  i.e., 


P{^Sg)  = 


Pi  if  0  <  z  < 
P2  \i  h  <  z 


In  the  upper  layer  (i.e.,  for  0  <  z  <  h),  the  pressure  field  can  be  written  in  terms 
of  the  sum  of  a  residue  series  and  a  branch  line  integral.  In  the  present  study,  the  branch 
line  integral  is  ignored  and  the  field  is  approximated  as  the  sum  of  the  terms  of  the  residue 
series.  The  explicit  form  of  the  series  is  given  as 


G(R„\R.,) 


sin  ctjZgg  sin  ajZ}gH^\rAjki  —  a 


2h 


7i  \  gjn2  Q/  ./itanctj/i 

h  a-i  a‘^-  J  •/ 


where 


b  =  ^,ki 

P2 


UJ 

Cl 

5 

7 

^2  “  — 5 

71=  — 

Cl 

^2 

C2 

and  where  the  terms 
dispersion  relation 


Xj  —  are  poles  of  G  which  are  located  at  values  satisfying  the 


tan(ah)  = 


ah 


by/h^kf(l  -  n2)  -  (ahy 


In  the  foregoing  equations,  the  subscript  i  refers  to  the  properties  of  the  upper  layer 
and  the  subscript  2  refers  to  the  properties  of  the  lower  layer.  The  term  ki  is  the  reference 
wave  number  mentioned  previously,  i.e.,  ki  —  ko.  Furthermore,  only  the  real  solutions  of  the 
dispersion  relation  are  to  be  considered  in  the  present  instance.  The  complex  solutions  will 
correspond  to  nonpropagating  modes  which  decay  exponentially. 

In  the  nomenclature  developed  in  the  foregoing  sections,  one  has  that 

(f)Az)  = - sin(Q;jz)  , 
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and 


C,  = 


-2 


hi  —  152^1(1; 

h  a 


sin^  Oj/i  tan  a^h 


The  coefficients  \j  and  Cj  have  been  computed  off-line  and  are  stored  in  a  disk  file  for 
use  with  the  Green’s  function  program. 

A.9  UNIFORM  WAVEGUIDE  WITH  HARD  BOTTOM 

In  the  case  of  this  environment,  a  surface  liquid  layer  of  depth  h  lies  over  a  hard 
bottom,  i.e.,  the  normal  derivative  of  the  Green’s  function  has  value  zero  at  the  bottom. 

In  the  liquid  layer,  the  pressure  held  can  be  written  in  terms  of  the  sum  of  a  residue 
series  (without  a  branch  line  integral).  In  the  present  study,  the  branch  line  integral  is 
ignored  and  the  held  is  approximated  as  the  sum  of  the  terms  of  the  residue  series.  The 
explicit  form  of  the  series  is  given  as 


G{Rjg\Rsg) 


2h 


E 

j=l 


sin  cXjZsg  sin  ajZfgllQ 


(2) 


[ry/W- 


where  the  terms  \j  =  aj  are  poles  of  G,  which  are  located  at  values  satisfying  the  dispersion 
relation 


cos{ah)  =  0  . 


The  complex  solutions  of  the  dispersion  relation  will  correspond  to  nonpropagating  modes 
that  decay  exponentially. 

In  the  nomenclature  developed  in  the  foregoing  sections,  one  has  that 

=  —  sin(Q;j2)  , 


and 


The  coefficients  Xj  and  Gj  have  been  computed  offline  and  are  stored  in  a  disk  hie  for  use 
with  the  Green’s  function  program. 
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The  foregoing  mode  shape  can  be  cast  into  a  form  that  is  particularly  useful  for 
interpreting  the  inter-modal  target  strengths  defined  previously.  Specifically,  note  that 


The  partial  source  strength  used  for  the  calculation  of  D  in  equation  (A.  12)  is  thus 


2  (  \  / ^j’f'srf  '^z  ~ 


where  use  has  been  made  of  the  relationship 


^srj  —  ^ srf  ■ 


This  result  is  equivalent  to  having  the  target  ensonified  by  two  free  field  waves,  one  of  which  is 
traveling  in  an  upward  direction  and  the  other  of  which  is  traveling  in  a  downward  direction. 
The  depression/elevation  angle,  0 depression /elevation  is  given  by  the  expression 


COs(^d depression /elevation^  —  •^1  Xj  j  /Jq 


A.IO  COORDINATE  SYSTEM  COORDINATION 

The  formulations  of  the  propagation  and  structural  problems  each  lend  themselves  to 
expression  in  a  particular  coordinate  system.  The  systems  are  cylindrical,  in  the  case  of 
the  propagation  problem,  with  the  origin  lying  at  the  surface  of  the  water  and  Cartesian, 
in  the  case  of  the  structural  problem,  with  the  origin  lying  somewhere  inside  of,  or  on,  the 
surface  of  the  shell.  The  cartesian  coordinates  to  which  reference  is  made  for  the  structural 
problem  are  called,  in  the  current  CHIEF/MARTS  AM  usage,  “global  coordinates.”  This 
distinguishes  them  from  the  “local  coordinates,”  which  are  used  in  the  individual  finite 
elements.  The  assumption  is  made  that  the  Xi  axis  lies  along  the  axis  of  symmetry  of  the 
body  being  considered.  The  X2  and  0:3  axes  are  chosen  such  that  Xz  points  in  the  vertical 
downward  direction,  x^  lies  in  the  horizontal  plane,  and  the  set  forms  a  right-handed  triad. 
In  the  case  of  the  cylindrical  system,  the  z  coordinate  (depth)  is  taken  to  be  a  translated 
version  of  xz  in  the  Cartesian  system.  The  radial  variable  in  the  cylindrical  system  is  given 
by  the  expression 


=  \Ja 


-f  X 


An  azimuthal  angle  variable  is  needed  to  specify  the  aspect  of  the  shell  relative  to 
the  field  point  location  (source  or  receiver).  The  convention  chosen  is  the  one  currently  used 
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by  the  CHIEF/MARTSAM  software  to  specify  target  aspect  relative  to  an  incident  plane 
wave.  A  second  angle  (elevation)  is  ignored  in  this  study  since  it  is  not  particularly  useful 
in  the  context  of  a  refractive  environment. 

Given  the  cartesian  coordinates  (xi, .7:2, 2:3)  of  a  point  on  the  wetted  surface  of  the 
shell,  then  the  cylindrical  coordinates  of  that  point  (relative  to  the  field  point)  are  given  as 


r  =  (xi  -  Tq  cos  4>a0y  +  {^2  “  ’’O  siu  (j^aoy 
z  =  zq  X3  ,  (A. 14) 


where  tq  is  the  cylindrical  range  from  the  field  point  to  the  origin  of  the  Cartesian  system, 
Zq  is  the  depth  (in  the  cylindrical  system)  of  the  origin  of  the  Cartesian  system,  and  </)ao 
is  the  aspect  of  the  target  relative  to  the  line  joining  the  field  point  and  the  Xi  axis  of  the 
cartesian  coordinate  system.  In  previous  sections  (A. 8  and  A. 9),  it  is  equation  (A. 14)  that 
is  used  to  calculate  coordinates  for  use  in  the  generation  of  coefficients  for  use  in  equation 
(A. 6)  and  equation  (A. 13). 

A.ll  APPROXIMATIONS  IN  SURFACE  PRESSURES  AND  VELOCITIES 

A  particular  approximation  is  used  in  the  calculation  of  the  surface  pressures  and 
velocities  that  is  justified  in  the  particular  cases  that  will  be  of  interest  in  the  present 
study.  The  residue  series  representation  of  the  propagation  Green’s  function  has  regions 
in  depth  and  range  for  which  it  does  not  converge.  While  these  regions  are  not  important 
from  the  standpoint  of  the  calculation  of  the  incident  pressure  from  a  distant  source  and 
the  far-held  scattered  pressures,  they  are  important  from  the  standpoint  of  the  solution  of 
the  integral  equation  dehning  the  surface  pressures  and  velocities  on  the  excited  shell.  It 
was  deemed  reasonable  to  use  the  free-held  propagation  Green’s  function  in  the  calculations 
involved  in  solving  the  integral  equation.  The  heuristic  justihcation  derives  from  the  notion 
that  the  distances  between  points  on  the  shell  are  much  smaller  than  their  distances  from 
corresponding  points  on  the  image  shell,  and  also  from  the  notion  that  refractive  effects 
at  these  distances  will  be  relatively  small.  This  approximation  can  also  be  cast  into  the 
form  of  a  low-order  term  in  a  proposed  sequence  of  approximations  to  the  actual  solution. 
This  sequence  is  developed  in  the  following  analysis.  Write  the  differential  equation  for  the 
Green’s  function  in  the  form 


CiG)  =  6-6_  , 


where  the  differential  operator  C  represents  the  wave  equation  in  a  refractive  environment 
that  includes  not  only  the  half-space  in  question,  but  also  for  an  image  half-space  above, 
i.e.,  the  composite  medium  has  a  sound  speed  profile  that  is  symmetric  in  depth  about  the 
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depth  origin  in  the  actual  medium.  In  addition,  the  terms  6  and  represent  sources  at 
the  actual  source  point  and  at  the  image  source  point  above  the  surface,  respectively.  The 
change  in  sign  for  the  image  point  presence  of  a  pressure  release  boundary  as  opposed  to  a 
rigid  boundary. 

The  foregoing  equation  can  be  rewritten  in  terms  of  a  wave  operator  for  a  uniform 
medium  and  a  residual  operator  embodying  the  depth-dependent  speed  of  sound.  The  result 
is 


(Co  +  6C)G  =  6  -  S.  , 


where  Co  is  a  differential  operator  for  an  approximate  uniform  medium,  and  6C  embodies 
the  refractive  effects  of  the  nonuniform  sound  speed  profile.  This  equation  is  written  in  a 
somewhat  more  convenient  form  as 


CoG=6-6_-  6CG  . 

Using  the  inverse  operator,  Co,  one  has  the  result  that 


GooCoG  =  Goo6-Goo6_-Goo  8CG , 


where  the  implied  integration  is  over  the  entire  domain  of  definition  of  G  and  C. 

From  the  definition  of  Go,  one  has  the  problem  cast  into  the  following  integral  equation: 


G  —  Go  —  Go—  —  Go  0  6CG.  (A. 15) 

The  problem  of  solving  for  the  Green’s  function  is  hence  cast  into  the  form  of  solving 
the  foregoing  problem.  Under  certain  circumstances,  the  problem  can  be  solved  by  the 
method  of  Picard  iteration.  Hence,  one  generates  a  sequence  of  approximations,  G[„],  to  the 
function  G  with  the  hope  that  the  sequence  converges  to  the  desired  function.  Specifically, 
one  starts  with  Go  and  Go-  and  one  constructs  the  sequence 


G[i]  —  Go  —  Go-  —  Go  o  SCGq 

G[2]  =  Go  —  Go-  —  Go  o  8CG[i\ 

G[n-\-\]  =  Go  —  Go-  —  Go  0  8CG[n] 
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This  sequence  of  approxinaations  to  the  Green’s  function  will,  under  certain  circumstances 
(unknown  at  this  point),  converge.  It  will,  for  no  particularly  good  reason,  be  assumed  to 
converge.  In  this  context,  the  approach  to  be  used  in  the  current  study  can  be  thought  of  as 
the  0-th  order  approximation.  While  the  higher  order  approximations  have  not  been  used 
in  this  study,  the  value  of  this  formulation  is  that  it  provides  a  constructive  approach  for 
calculating  each  term  of  the  sequence. 

A  variant  of  the  lowest  order  approximation  which  is  of  particular  interest  is 

G[o'](jRo6s|-Rsp)  —  Go  —  Go--  (A. 16) 

This  approximation  will  capture  important  multiple  scattering  effects  for  a  target  near  the 
pressure  release  boundary. 

A.12  APPROXIMATE  MULTIPLE  SCATTERING  GREEN’S  FUNCTION  FOR 
SURFACE  INTEGRAL  EQUATION 


The  Surface  Integral  Equation  calls  for  the  use  of  the  full  refractive  Green’s  function 
for  use  in  the  kernel.  Difficulties  with  a  convenient  representation  suggest  the  use  of  simpler, 
approximate  forms.  A  method  for  partially  implementing  the  approximate  scheme  in  the 
previous  section  is  detailed  here. 


Let  the  vectors  Rfp  and  Rsp  be  the  field  and  source  points  for  the  Green’s  function  in 
the  following  free-space  boundary  value  problem 


+  kl)Go{Rfp\Rsp)  —  0  in  U  , 


and 


lim  Vr 


I  dGo 

dr 


+  ikoGol  —  0  , 


where  Vp  =  \Rfp  -  R^p],  ko  =  uj/cq- 


In  the  present  example,  let  Rfp  =  Rfg  =  Rf.  We  then  have  the  result 


GoiRf\Rsp){Vi^  +  k\zf))G{Rf\R,g)  -  G{Rf\Rsg){Vj,^  +  k^,)Go{RfAp)  =  0  , 


or 


(kHz,)  -  e„)Go{R,\R„)G(R,\R„)  =  0  .  (A.17) 

The  following  identities  are  useful  in  simplifying  the  foregoing  expression: 

Vfi,  •  [GoiRf\Rsp)yR,G{Rf\Rsg)]  =  Vn^GoiRfAp)  ■  VR,G{Rf\R,g)  + 

Go{Rf\Rsp)^i^G{Rf\Rsg)  ^  0  , 
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and 


Rj  '  [G{Rf\Rsg)'^ RfGo{Rf\Rsp)]  —  ^RfG{Rf\Rsg)  • '^RjGo{Rf\Rsp)  + 

G{Rf\Rs,)Wi^Go{Rf\Rsp)  =  0  . 

Hence,  we  can  rewrite  equation  (A.  17)  in  a  simplified  form  as 

VR^-[GoiRf\Rsp)yR,GiRf\R,,)  -  G{Rj\R,g)VR^GoiRf\Rsp)]  + 

{k^{zf)  -  kl)Go{Rf\Rsp)G{Rf\R,,)  =  0. 


This  expression  can  be  integrated  over  the  volume  that  excludes  the  shell  and  its 
interior,  a  tiny  sphere  of  radius  e  centered  at  i?sp,  and  another  tiny  sphere  of  radius  e 
centered  at  R^g.  Since  there  are  no  sources  in  this  volume,  one  has  the  result 

/  {Vh,  •  [Go{Rf\Rsp)VRfGiRj\Rsg)- 

G{R}\R,g)VR,Go{Rj\Kp)\  +  [k\zi)  -  kl]Go{RARsp)G{Ri\Rsg)}dV  =  0. 

This  is  readily  converted  into  a  surface  integral  of  the  form 

[  [Go{Rf\Rsp)^RjGiRf\R,g)  -  G{Rf\R,g)VnjGo{Rf\Rsp)]  ■  n^dS 

J  S^g‘^Scp~\-SuPPBR’\-SLOWER 

+  /  [k\zf)  -  kl]Go{R}\Rsp)G{Ri\Rsg)dV  =  0  ,  (A.18) 

JV  —  Vcg  —  Vep 

where  Ufi  is  the  outward  unit  normal  to  the  fluid. 


The  integrals  over  the  surfaces  S(,g  and  S^p  have  particularly  simple  limiting  forms. 
Note  that 

I  GoiRf\Rsp)^R,G{Rj\R,g)-nfidS  =  Go{Rsg\Rsp)  [  VR^G{Rf\R,g)  ■  HfidS 
^  S^g  J  S^g 

-  -  r  d  /e-**'' 


=  Go{Rsg\Rsp)  -(47re^) 


de  \  47re 


—  Go{Rsg\Rsp) 
^  Go{Rsg\Rsp) 


—  (47re^)  f  —ik- 


47r€  47re^ 


Similarly,  one  has 

f  G{Rf\Rsg)'V R^Go{Rf\Rsp)  ’  TifidS 

'd  Sc-n 


G{Rsp\Rsg)  j  RjGo{Rj\Rsp)  '  fifldS 
"  Srn 


d 


-ike 


—  G{Rsp\Rsg) 
=  G{Rsp\Rsg)  • 


— (47re^)  —ik 


^-ike  ' 


47r6  47re^ 
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Note  also  that 


/  G{Rj\Rsg))\7R^GoiRf\Rsp)-nf,dS  ^  \^RjGo{RJRsp)\  [  GiRf\Rss)e-nj,dS 

JScp 

—  |Vij^G'o(^s5l-Rsp)|T —  /  e-njidS 

aTT  €  */  Sirr> 


0. 


Similarly,  one  has 

I  Go{Rf\Rsp))'^ RfG{Rf\Rsg)  ■  UftdS 

Sen 


\VR^G{Rsp\Rsg)\  j  Go{Rj\Rsp)e-npdS 

S  eg 

\V  RjG{Rsp\Rsg)\- —  /  e-UfidS 

“lTT  C  «/  Sen 


0. 


These  four  approximations  are  exact  in  the  limit  as  e  shrinks  to  0.  The  surface  integral 
can  thus  be  written  as 


Go{Rsg\Rsp)  “  G{Rsp\Rsg)-{- 

/  [Go{RuPPER\Rsp)'^RupPER^i^UPPER\Rsg)  “ 

Supper 

G(i?[/pp£'/j|i255)Vfiypp^^Go(-Rc/pp£;H|-R5p)]  •  nftdS  + 
r  _  _ 

[Go{RLOWER\Rsp)'^RlowER^i^LOWER\Rsg)  — 


•^Slower 

G{RLOWER\Rsg)'^  R, 

The  surface  integral  can  thus  be  rewritten  as 
Go{Rsg\Rsp)  =  G{Rsp\Rsg)  — 


^Go{Rlower\Rsp)]  ■  njidS  =  0. 


/  [Go{RuPPER\Rsp)'^RuPFER^i^UPPER\Rsg)  — 

^Supper 

G{RupPER\Rsg)^RuppBRGo{RuPPER\Rsp)\  '  PpdS 
-  /  [Go{RlOWEr\Rsp)^ RLOWERG{RLOWER\Rsg)  “ 

^Slower 

ER\^sgy^ Rlow er^^^^LOW Er\^s'p)\  *  RfldS. 


(A.19) 


(A.20) 
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A  change  of  notation  can  be  invoked,  and  the  symmetry  of  the  Green’s  function  can 
be  used  to  cast  the  problem  into  a  form  familiar  to  those  working  in  scattering  theory 

Go{^Robs\^sp)  — 

/  [Go{RuPPER\Rsp)^RuppER^i^obs\RuPPER)  — 

•^Supper 

— ^  ^  , 

G{Robs\RuPPER)'^RuppERGo{RuPPER\Rsp)]  '  RfidS 

—  /  [Go{RloWEr\Rsp)^  RLowERG{Robs\RLOWER)  — 

'^Slower 

G{Robs\RLOWER)'^RLOWERGo{RLOWER\Rsp)]  '  flfidS 
+  l^{k\zf)  -  Po)Go{RfAp)G{Rf\Rsy)dV  , 

where  Robs  is  substituted  for  the  term  Rsg^ 

Using  the  pressure  release  upper  boundary  condition  and  the  rigid  lower  boundary 
condition  allows  one  to  cast  the  equation  into  the  following  form: 


Go(^Robs\Rsp^  ' —  ^(^Robs\Rsp} 

/  [Go{RuPPER\Rsp)^RupPER^i^obs\RuPPER)]  *  flfidS 

'^Supper 


+  /  \S RLowERGo{RLOWER\Rsp)G{Robs\RLOWER)]  '  RjidS 

•^Slower 

+  f  {k^zj)  -  kl)GoiRf\Rsp)G{Rf\R,g)dV 

J  V 


or 


G{Robs\Rsp)  —  Go{Robs\Rsp)~\~ 

/  [Go{RuPPER\Rsp)'^RuppER^i^obs\RuPPER)]  '  RfldS 

J Supper 

~  /  RLOWER^^iRi^OWER\Rsp)G{Robs\RLOWER)]  ’flfidS 

^Slower 


-  j^(e{z,)  -  kl)Ga(Rt\R.p)G(R,\R.,)dV  . 


(A.21) 


Equation  (A.21)  is  equivalent  to  equation  (A. 15).  The  interpretation  of  the  above 
result  is  that  the  resulting  field  consists  of  the  direct  free-space  arrival  from  the  source 
represented  by  Go 5  ^  term  scattered  off  the  pressure  release  surface  and  represented  by  the 
surface  integrals,  and  a  volume  term  due  to  refraction. 
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Consider  the  special  case  of  a  uniform  environment  with  a  pressure  release  upper 
surface  and  a  rigid  lower  surface.  The  foregoing  equation  thus  becomes 

/  [Go{RuPPER\Rsp)'^RfiG{Robs\RuPPER)]  '  nyppERdS 

•^Supper 


+  /  [^RfiGo{RLOWER\R  sp  )G{Robs\RLOWER)]  ■  R-LOWERdS. 

'^Slower 


This  equation  can  also  be  written  in  terms  of  image  sources  as  follows 

Go{Robs\Rsp)  =  G{Robs\R  sp)  ^0 (-^o6s  I -^sp  upper  image)  ^0 (-^o6s  |-^sp  lower  image)  5  (^A..22) 


where  R 
surface. 


sp  image 


is  the  location  of  the  image  of  the  source  point  above  the  pressure  release 


A  special  set  of  coordinates  will  be  chosen  with  the  origin  at  the  phase  center  of  the 
target.  The  pressure  release  surface  will  have  a  z- coordinate  of  zipper  and  the  lower  rigid 
boundary  will  have  a  z-coordinate  of  zpowER-  Using  equation  (A. 22),  the  Green’s  function 
(equation  A.21)  is  hence  written  approximately  as  follows 


^?[o^]  (/?o6s  l^sp)  — 


0  Es 


R,^ 


•sp  upper  image 


47r|i?o6s  -  R. 


spl 


47r|^, 
e 


+ 


obs  ^sp  upper  image  \ 
—ik\Robs  —  E 

sp  lower  imagel 

irlR^bs  -  R  sp  lower  image  \ 


(A.23) 


where 


\R 


obs 


R 


spl 


i^^obs  ^sp )  T  (^0^  Vsp)  T~(^o6s  ^sp)  5 


j^o6s 


sp  upper  image 


and 


l^obs 


sp  lower  imagel 


—  \/[^obs  ~  XgpY  +  {Vobs  —  PspY  +  i^obs  ”  {^WPPER  —  5 


It  is  equation  (A.23),  which  is  an  implementation  of  equation  (A. 16),  that  is  used  in  place 
of  Go  in  the  Surface  Integral  Equation  (A. 6). 

An  alternative,  which  is  valid  only  for  the  approximation  G[o'](.Ro6s|.Rsp),  is  to  use  an 
image  target  with  the  pressure  and  velocities  being  of  opposite  sign  to  those  of  the  true 
target.  Equation  (A. 6)  can  be  manipulated  as  follows  to  yield  a  form  that  is  particularly 
easy  to  implement: 

-P(Rsg\Rsp)  =  G{Rsp\Rsg)-\- 

f  [P{Rsh\Rsp)'^ R,hG[0']iRsh\Rsg)  -  G[0']{Rsh\Rsg)'^ R,hP{Rsh\Rsp)]  ■  RshdS  , 

S  oh 
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or 


^P(R.,\R„)  =  G4R„\R„)  -  G.{R 

sp  image  \R 

sg  )+ 

/  [P{Rsh\Rsp)'^ R,h^oiRsh\Rsg)  -  Go{Rsh\Rsg)'^ R.hP{Psh\Rsp)]  '  nshdS  — 
f  [PiRshAp)^RsRGo{Rsk\Rsg)  -  Go{Rsh\Rsg)^R.RP{Rsh\Rsp)]  ■  ftshdS  , 

'^^sh  image 

where  the  functions  G+{Rsp\Rsg)  and  G-{Rsp  image\Rsg)  are  the  source  terms  (Green’s  func¬ 
tions)  for  the  environments  which  have  the  sound  speed  profile  reflected  up  above  the  pressure 
release  surface,  i.e.,  a  mirror  image  environment  lying  above  the  surface. 

The  specific  computer  implementation  would  be  as  follows.  Solve  two  separate  prob¬ 
lems  and  add  the  solutions: 

P{Rsg\Rsp)  =  P+{Rsg\Rsp)  -  P-iRsg\Rsp)  , 


where 


-P+{Rsg\Rsp)  —  G+{Rsp\Rsg)  + 


and 


/  [P+{Rsh\Rsp)'^RshGo{Rsh\Rsg)  —  Go{Rsh\Rsg)^ Rshd^+i^sh\Rsp)]  '  nghdS  — 
JSsh 

/  [P+{Rsh\Rsp)'^ RshGo{Rsh\Rsg)  —  Go{Rsh\Rsg)'’^ R,hd^+{^sh\Rsp)]  '  RshdS 

^sh  image 


^P—(,Rsg\Rsp')  —  G  —{^Rsp  image\Rsg^P 

f  [P-{Rsh\Rsp)'^ R,hGoiRsh\Rsg)  —  Go{Rsh\Rsg)'^ RshP-iPsh\Rsp)]  '  dshdS  — 

I  [P^{Rsh\Rsp)'^ R^hGo{Rsh\Rsg)  —  Go{Rsh\Rsg)'^ R.hd^-i^sh\Rsp)]  '  RghdS  . 

JSsh  image 


A  separate  target  must  be  placed  at  the  target  image  location,  and  the  appropriate 
change  in  sign  must  be  made  in  its  contributions  to  the  integral  equation.  It  may  require 
reprogramming  the  integral  equation  solver  in  CHIEF  to  put  in  this  change  of  sign.  Since 
the  “imaging”  effect  is  only  due  to  the  particular  representation  used  for  the  propagation 
Green’s  function,  no  software  changes  should  be  needed  with  respect  to  the  usage  of  the 
elastic  target  effects.  Note  that  the  driving  functions  for  the  equations,  G-i.{Rsp\Rsg)  and 
G-{Rsp  image\Rsg)-,  ^re  Gieeii’s  functioHS  for  the  reflected  environment,  rather  than  for  the 
original  half  space;  thus,  care  must  be  taken  to  use  the  original  modes  correctly  when  the 
target  and  source  are  not  in  the  same  half-space. 
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From  symmetry  considerations,  one  has  that 

P^(^Rsg\Rsp^  P—^Rsg  image\Rsp)  • 

This  means  that  one  need  only  solve  one  of  the  two  foregoing  integral  equations  and  then 
correctly  combine  the  solution  for  the  true  target  with  the  solution  at  the  appropriate  image 
points  of  the  image  target  in  order  to  get  the  full  solution. 

A.  13  SPECIAL  PARTIAL  SUMS  OF  INTEREST 

Equation  (A. 22)  can  be  rewritten,  for  the  case  wherein  the  effect  of  the  upper  surface 
is  to  be  neglected,  as  follows: 


G(^Robs\Psp)  —  Goi^RobslRsp) 

—  /  RfiGo{RLOWER\Rsp)G{Robs\RLOWER)]  •  nLoWERdS. 

'’Slower 

Note  that  the  equation  has  now  been  cast  into  the  form  of  a  fixed  point  problem. 
In  the  event  that  the  right-hand  side  is  a  contraction  mapping  and  that  one  can  locate 
a  starting  point  in  the  domain  of  attraction  of  the  fixed  point,  then  a  candidate  solution 
technique  is  Picard  iteration.  In  particular,  the  problem  then  becomes  equivalent  to  the 
following  sequence  of  problems: 


G^^\Robs\Rsp)  =  Go{Roi,s\Rsp) 

~  /  ['^RfiGo{RLOWER\Rsp)Go{Robs\RLOWER)]  '  nioWERdS  , 

St. nw-RTf 


and 


G^'^\Robs\Rsp)  =  Go{Robs\K) 

~  L  [^RfiGo{RLOWER\Rsp)G^^\Robs\RLOWER)]  •  niOWERdS. 

•d  Slower 

The  n-th  iteration  is  given  by 

=  Go{Robs\Rsp) 

—  /  [^RjiGo{RlOWEr\R  sp  )G^^-^\R 

obs  \Rlower)]  ■  niowERdS. 

^Slower 

Use  is  to  be  made  of  equation  (A. 22),  which  can  be  rewritten  (neglecting  the  upper  surface) 
as 

GiRobs\Rsp)  =  Go{Rohs\R  sp)  T  G()(.RQbs\Rsp  lower  imaged  • 
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The  foregoing  equations  yield  the  result 

I )  * —  ^o(-^o65  |-^sp) 

~  IsiowE^^^fi^^^^^OWERlRsp)  {Go{Robs\RLOWER)  + 
Go{Robs  image  \Rlower)\]  •  nwwERdS  , 
so  the  first  iterate  is  given  as 

g^HR  obs  l-Rsp)  —  Go(Robs\Rsp) 

—  /  [^RfiGo{RLOWER\Rsp)Go{Robs\RLOWER)]  '  nioWERdS 

•'Slower 

=  Go(-Ro6s|-R  sp) 

\  f  — *  >  >  > 

—  X  /  ['^Rj,Go{RLOWER\Rsp)G{Robs\RLOWER)]  ‘ 'dlOWERdS 

^  ’d  S T.nWK  R 


1 


Go(^Robs\Rsp)  T  2^^Q{Robs  lower  image\Rs'p)  • 

The  second  iterate  is  given  by 
obs  l-^sp)  —  Go{Robs\Rsp) 

—  /  RfiGo{RLOWER\R  sp  )G^^\R 

obs  \Rlower)]  •  ntowERdS 

^Slower 

~  Go{Robs\R  sp) 

—  f  RpGo{RLOWER\Rsp)Go{Robs\R  sp)  H”  q^o(-^o6s  lower  image  \R  5p)]  ’  'l^LOWER^'S 

d  Slower  ^ 

—  Go{Rohs\R  sp  )  +  -^GoiRobs  lower  image\Rsp  )  +  -^Go{Robs  lower  image\Rsp)  ■ 

The  n-th  iterate  has  the  simple  form 

G^  ^ (-Ro6s  |-^sp)  \  ^  V  7^  I  ^O^Robs  lower  image\Rsp) 

[j=0^  J 


or 


f  ^  1  1 

o65i-^5p)  A  y  q-  J  GQ(^RQbs\Rsp  lower  imaged  • 

lj=0^  J 


Alternatively,  equation  (A. 22)  can  be  rewritten,  for  the  case  wherein  the  effect  of  the 
lower  surface  is  to  be  neglected,  as  follows: 


G{Robs\Rsp)  —  Go{Robs\Rsp)~\- 

/  [Go(Rupper\R  sp  )^Rf,G{Robs\RuPPER)]  •  duppEpdS  . 

■J  Supper 
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Note  that  the  equation  has  now  been  cast  into  the  form  of  a  fixed  point  problem.  In 
the  event  that  the  right-hand  side  is  a  contraction  mapping,  and  one  can  locate  a  starting 
point  in  the  domain  of  attraction  of  the  fixed  point,  then  a  candidate  solution  technique 
is  Picard  iteration.  In  particular,  the  problem  then  becomes  equivalent  to  the  following 
sequence  of  problems: 


G^^KRobslRsp)  =  GoiRo,sAp)+ 

/  [Go{RuPPER\Rsp)'’^RjiGo{Robs\RuPPER)]  '  ^UPPERdS  , 

^^UPPER 


and 


G^‘^\Rohs\Rsp)  =  Go{Robs\Rsp)+ 

/  {Go{RuPPEr\Rsp)^ RjtG^^\Robs\RuPPER)]  '  RUPPERdS 

J  StT  IDTDTP  -D 


The  n-th  iteration  is  given  by 

g^^\r 

ohs  \Rsp)  =  Go{R  obs  |^5p)  + 

/  [Go{RuPPER\Rspy^RfiG^'^~^\Robs\RuPPER)]  '  RUPPERdS 

Supper 


Use  is  to  be  made  of  equation  (A. 22),  which  can  be  rewritten  (neglecting  the  upper 
surface)  as 

Gr{Robs\Rsp)  =  Go{Robs\R  5p)  G()(^RqI)s\Rsp  upper  image')  * 

The  foregoing  equations  yield  the  result 

G{Robs\Rsp)  =  Go(Robs\Rsp) 

+  /supp£H  Go ( Rupper  I Rsp )  { Go ( .Ro6s  I Rupper)  - 

GoiRobs  image  \RuPPEr)^]  ■  RVPPERdS  , 


so  the  first  iterate  is  given  as 

G^^\R,bs\Rsp)  =  Go(^.;,.|^.p) 

+  /  RfiGo{RuPPER\R  sp  )Go{Robs\RuPPER)]  •  RUPPERdS 

Supper 

=  Go{Robs\R 

1  /*  _  -*  -* 

+  9  /  [^RfiGo{RupPER\R  sp  )G{Robs\RuPPER)]  •  nuPPERdS 

Z  SrjPPF.R 


1 


=  Go{Robs\R  sp  )  2^Gq(^RqI)^  upper  image\Rsp)  * 
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The  second  iterate  is  given  by 


obs  l-Rsp)  —  Go{Robs\Rsp) 

+  [  [VRj,Go{RuPPER\Rsp)G^^\Robs\RuPPER)]  ■  nUPPERdS 

^Supper 

=  Go{Robs\R  5p) 

+  /  ['^RfiGo{RuPPER\Rsp)Go{Robs\R  sp)  Gq(^Ro})S 

upper  image\Rsp  )]  •  njjppERdS 

-'Supper  ^ 

=  Go(^o65|^  sp)  2  Goi^Robs  upper  image\Rsp  )  ^G{}iyRobs  upper  image\Rsp^  • 

The  n-th  iterate  has  the  simple  form 
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This  final  equation  ends  the  supporting  analysis  for  the  main  body  of  the  report. 
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